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ABSTRACT 

This paper deals with the free torsional vibrations of doubly symmetric thin-walled beams of open section and 
resting on Winkler-Pasternak continuous foundation. A general dynamic stiffness matrix is developed in this paper which 
includes the effects of warping and Winkler-Pasternak foundation on the torsional natural frequencies. The resulting highly 
transcendental frequency equations for a simply supported and guided-end conditions are solved for varying values of 
warping Winkler and Pasternak foundation parameters on its frequencies of vibration. A new MATLAB code was 
developed based on modified BISECTION method to solve the highly transcendental frequency equations and to 
accurately determine the torsional natural frequencies for various boundary conditions. Numerical results for natural 
frequencies for various values of warping and Winkler and Pasternak foundation parameters are obtained and presented in 
graphical form showing their parametric influence clearly. 
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INTRODUCTION 

The problem of a beam on an elastic foundation is important in both the civil and mechanical engineering fields, 
since it constitutes a practical idealization for many problems (e.g. the footing foundation supporting group of columns (as 
shown in Figure 1) etc. The concept of beams and slabs on elastic foundations has been extensively used by geotechnical, 
pavement and railroad engineers for foundation design and analysis. 
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The analysis of structures resting on elastic foundations is usually based on a relatively simple model of the 
foundation’s response to applied loads. A simple representation of elastic foundation was introduced by Winkler(as shown 
in Figure 2) in 1867 [1]. The Winkler model (one parameter model), which has been originally developed for the analysis 
of railroad tracks, is very simple but does not accurately represents the characteristics of many practical foundations. 




Figure 2: Deflections of Winkler Foundation under Uniform Pressure Q 

In order to eliminate the deficiency of Winkler model, improved theories have been introduced on refinement of 
Winkler’s model, by visualizing various types of interconnections such as shear layers and beams along the Winkler 
springs [1] (Filonenko -Borodich (1940) [2]; Hetenyi (1946) [3]; Pasternak (1954) [4]; Kerr (1964) [5]). These theories 
have been attempted to find an applicable and simple model of representation of foundation medium. To overcome the 
Winkler model shortcomings improved versions [6] [7] have been developed. A shear layer is introduced in the Winkler 
foundation and the spring constants above and below this layer is assumed to be different as per this formulation. The 
following figure shows the physical representation of the Winkler-Pasternak model. 



MpjlMpPf 

Winkler * 

springs Shearing 

layer {jmsfcniak springs) 



Figure 3: Winkler-Pasternak Model 

The vibrations of continuously- supported finite and infinite beams on elastic foundation has wide applications in 
the design of aircraft structures, base frames for rotating machinery, railroad tracks, etc. Quite a good amount of literature 
exists on this topic, and valuable practical methods for the analysis of beams on elastic foundation have been suggested. [8- 
11]. Karnes wara Rao et al [12-14] studied the problem of torsional vibration of long, thin- walled beams of open sections 
resting on Winkler-type elastic foundations using exact, finite element and approximate expressions for torsional frequency 
of a thin- walled beams and subjected to a time-invariant axial compressive force. 

It is well known that a dynamic stiffness matrix is mostly formed by frequency-dependent shape functions which 
are exact solutions of the governing differential equations. It overcomes the discretization errors and is capable of 
predicting an infinite number of natural modes by means of a finite number of degrees of freedom. This method has been 
applied successfully for many dynamic problems including natural vibration. A general dynamic- stiffness matrix of a 
Timoshenko beam for transverse vibrations was derived including the effects of rotary inertia of the mass, shear distortion, 
structural damping, axial force, elastic spring and dashpot foundation [15]. Analytical expressions were derived for the 
coupled bending-torsional dynamic stiffness matrix terms of an axially loaded uniform Timoshenko beam element [16-20] 
and also a dynamic stiffness matrix is derived based on Bernoulli-Euler beam theory for determining natural frequencies 
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and mode shapes of the coupled bending-torsion vibration of axially loaded thin-walled beams with mono-symmetrical 
cross sections, by using a general solution of the governing differential equations of motion including the effect of warping 
stiffness and axial force [21] and [22]. Using the technical computing program Mathematica, a new dynamic stiffness 
matrix was derived based on the power series method for the spatially coupled free vibration analysis of thin- walled curved 
beam with non-symmetric cross-section on Winkler and also Pasternak types of elastic foundation [23] and [24]. The free 
vibration frequencies of a beam were also derived with flexible ends resting on Pasternak soil, in the presence of a 
concentrated mass at an arbitrary intermediate abscissa [25]. The static and dynamic behaviors of tapered beams were 
studied using the differential quadrature method (DQM) [26] and also a finite element procedure was developed for 
analyzing the flexural vibrations of a uniform Timoshenko beam-column on a two-parameter elastic foundation [27]. 

Though many interesting studies are reported in the literature [8-27], the case of doubly- symmetric thin- walled 
open section beams resting on Winkler-Pasternak foundation is not dealt sufficiently in the available literature to the best 
of the author’s knowledge. 

In view of the above, the present paper deals with the free torsional vibrations of doubly symmetric thin-walled 
beams of open section and resting on Winkler-Pasternak continuous foundation. A general dynamic stiffness matrix is 
developed in this paper which includes the effects of warping and Winkler-Pasternak foundation on the torsional natural 
frequencies. The resulting highly transcendental frequency equations for a simply supported and guided-end conditions are 
solved for varying values of warping Winkler and Pasternak foundation parameters on its frequencies of vibration. A new 
MATLAB code was developed based on modified BISECTION method to solve the highly transcendental frequency 
equations and to accurately determine the torsional natural frequencies for various boundary conditions. Numerical results 
for natural frequencies for various values of warping and Winkler and Pasternak foundation parameters are obtained and 
presented in graphical form showing their parametric influence clearly. 

FORMULATION AND ANALYSIS 

Consider a long doubly- symmetric thin-walled beam of open-section of length L and resting on a Winkler-Pasternak type 
continuous foundation of Pasternak layer stiffness (/6 p ), Winkler torsional stiffness ( K w ) and undergoing free torsional 
vibrations. The corresponding differential equation of motion can be written as: 



d 4 0 , , d 2 0 

‘w ^^4 — (g C s + — P*p 




( 1 ) 



For free torsional vibrations, the angle of twist 0(z, t) can be expressed in the form, 

&(z, t) = x(z}e l{x}t 

In which x(z) is the modal shape function corresponding to each beam torsional natural frequency a). 
The expression for x(z) which satisfies Eq. (1) can be written as 

x(z) = A cos pL 4- B sin pL 4- C cosh aL 4- D sinh aL 



( 2 ) 



( 3 ) 
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In which aL and (3L are the positive, real quantities given by 



pL, aL = 



\+(K 2 +Yp 2 ) + V(tf 2 + r P 2 ) 2 + 4(A 2 ->v 2 ) 



(4) 




From Eq. (4) we have the following relation between aL and fLL 
(pL ) 2 - (aL ) 2 + tf 2 -Y P 2 

Knowing a and /?, the frequency parameter X can be evaluated using the relation 



l 2 - (aL) (pL) + y w 2 

(7) 

The four arbitrary constants A, B, C and D in Eq. (3) can be determined from the boundary conditions of the 
beam. For any single-span beam, there will be two boundary conditions at each end and these four conditions then 
determine the corresponding frequency and mode shape expressions. 

DYNAMIC STIFFNESS MATRIX 

In order to proceed further, we must first introduce the following nomenclature. The flange bending moment and 
the total twisting moment are given by M(z) andT(z). The variation of angle of twist 0 with respect to z is denoted 
by 6 (z) . Considering anti-clockwise rotations and moments to be negative, we have 



6(z) = d®fdz r hM(z) = — EC w d 2 0/dz 3 

, .. d 2 (p dtp 

t ™ = - ec "J? + gc *k 

Where EC W is termed as the warping rigidity of the section, 



* ( 10 ) 

Consider a uniform thin- walled I-beam element of length L as shown in the Figure 3. By combining the Eq. (3) 
and Eq. (8), the end displacements, 0(O)and 0(0) and end forces /iM( 0) and 7(0), of the beam atz = 0, can be expressed 
as 



( 8 ) 

(9) 



Index Copernicus Value: 3.0 - Articles can be sent to editor @impactjournals.us 




The Influence of Warping and Winkler-Pasternak Soil on the Torsional Vibrations 19 

of Thin- Walled Open Section Beams with Guided-End Conditions 



0(0) ^ 


1 


0 


1 !. 


J 0 


0(0) _ 


0 


a 


0 


□ p 


JiAf(O) [ 


EC w a 2 


0 


EC w fi 2 


0 


r(o) / 


- 0 


EC w afi 2 


0 


-EC w a 2 fl. 




Eq. (3.28) can be abbreviated as follows 



do 



s(0) = 1(0)1/ (12) 

In a similar manner, the end displacements 0(L),0(L) and end forces /iM(L), 7(L) of the beam atz = L, can be 
expressed as 

5(_L) = S(L)U 

Where 

f 0(i) \ 

sa) = 0(I) 

WiM 

V T(_L) ) 

[uY = {A B c D } 



(14) 

(15) 



[£(!)] 



c s | \ [j C □ U S 

—as ac jS S fiC 

EC w a 2 c EC w a 2 s -EC w fl 2 C -EC W {3 2 S 
-EC w af3 2 s EC w ap 2 c -EC w a 2 fiS -EC w a 2 pC. 



In which 

c = cos j BL,s = sin flL r C = cosh aL r S = sinh aL 



i 




Figure 4: Differential Element of Thin Wall I Beam 



(16) 



(17) 
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The equation relating the end forces and displacements can be written as 
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(18) 



By eliminating the integration constant vector U and designating the left end element as I and the right end as j, 
the final equation relating the end forces and displacements can be written as 



Ti/EC w 1 
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Eq. (19) is symbolically written as 



(19) 



(U = [/]{£/} 



( 20 ) 



In the Eq. (20) the matrix [/] is the ‘exact’ element dynamic stiffness matrix (DSM), which is also a square matrix. 
The elements of [/]are 



7n = H(a 2 + jff 2 }(cr£c + jffCs) 

7i2 = -H[{cc 2 - jfi 2 )(l - Cc) + 2ccfiSs] 
7i 3 = -H(a 2 + p 2 )(as + pS) 

7 W = ~H{a 2 + j 8 Z XC - c) 
hi = -Qi/apXa 2 + p 2 )(aSc - pcs) 
7 24 = (H/ap){a 2 + p 2 )(aS - fis) 



~ Jl4> J 33 — Jll> hi ~ JlZ> J 44 — JzZ 



H = afi/[2ap(l - Cc + ( jff 2 - cr 2 )£s] 



( 21 ) 



Using the element dynamic stiffness matrix defined by Eq. (20), one can easily set up the general equilibrium 
equations for multi-span thin-walled beams, adopting the usual finite element assembly methods. Introducing the boundary 
conditions, the final set of equations can be solved for eigenvalues by setting up the determinant of their matrix to zero 



METHOD OF SOLUTION 



Denoting the modified dynamic stiffness matrix as [J], we state that 

det\j\ = 0 



( 22 ) 
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The above equation yields the frequency equation of continuous thin- walled beams in torsion resting on Winkler- 
Pasternak type foundation. It can be noted that above equation is highly transcendental, the roots of equation can, therefore, 
be obtained by applying the bisection method using MATLAB code on a high-speed digital computer. 

A new MATLAB code was developed based on bisection method and was published in MATLAB Central official 
online library [29] which was cited and referred by few researchers. 

• These are some of the key highlights of the new MATLAB code 

• Primarily, it solves almost any given linear, non-linear & highly-transcendental equations. 

• Additional key highlight of this code is, the equation whose roots are to be found, can be defined separately in an 

m” file, which facilitates to solve multi-variable (example* 2 + y 2 + sin x + cos y = 0) highly-transcendental 
equations of any size, where the existing MATLAB codes fail to compute. 

• This code is made robust in such a manner that it can automatically save and write the detailed informationsuch as 
no. of iterations; the corresponding values of the variables and the computing time are automatically saved, into a 
(.txt) file format, in a systematic tabular form. 

• This code has been tested on MATLAB 7.14 (R2012a) [30] for all possible types of equations and proved to be 
accurate. 

Exact values of the frequency parameter X for various boundary conditions of thin- walled open section beam are 
obtained and the results are presented graphical form in this paper for varying values of warping, Winkler foundation and 
Pasternak foundation parameters. 

RESULTS AND DISCUSSIONS 



The approach developed in this paper can be applied to the calculation of natural torsional frequencies and mode 
shapes of multi-span doubly symmetric thin-walled beams of open section such as beams of I-section. Beams with non- 
uniform cross-sections also can be handled very easily as the present approach is almost similar to the finite element 
method of analysis but with exact displacement shape functions. All classical and non-classical (elastic restraints) boundary 
conditions can be incorporated in the present model without any difficulty. 

In the following, the end conditions for this problem can be written as 



0 = 0] 
M = Q> 
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e = o 


J fixed end 
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Considering a one element solution and applying these boundary conditions to Eq. (19) gives, 
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Table 1: Comparison of Obtained Frequency Equations with Ref. [32] 



s 

No. 


Boundary 

Condition 


Frequency Equations (Present Work) 


Frequency Equations 
(Ref. [32]) 


1 


Guided-guided 


sin(/?L) * sinh(aL) = 0 


sin(y„) = 0 


2 


Guided-free 


a 3 tan (/?L) + /? 3 tanh(aL) = 0 


tan (y n ) + tanh(y„. ) = 0 


3 


Guided-pinned 


2a 2 /? 2 


cosh(y n ) cos(y n ) = -1 


(a 4 4- /? 4 ) C0S ^ aL J cos ypLj t i - u 


4 


Guided-clamped 


/ 3 tan (/?L) -1- a tanh(aL) = 0 


tan(y„) + tanh(y n ) = 0 



The transcendental frequency equations (Eq. (4a), Eq. (4b), Eq. (4e) and Eq. (4f)) of Ref. [32], for generally 
restrained beams are observed to be same as equations obtained in the present work (see Table 1) of the present paper, but 
for only difference that the roots a and p are considered to be equal and have same sign and defined as (y n ). 

Furthermore, the equations for all the BC’s are also solved for values of warping parameter 0 < K < 20 and for 
various values of Winkler foundation parameter 0 < y w < 15 and values of Pasternak foundation parameter 0 < y p < 1 
and are presented in graphical form for the first mode. 



plot for infleunce of Winkler's foundation parameter on frequencey parameter A for 
guided-clamped beam 




plot for infleunce of winkler's foundation parameter on frequencey parameter Afor 
guided-free beam 
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plot for infleunce of Winkler's foundation parameter on frequencey parameter A for 
guided-guided beam 




pasternak foundation parameter 



plot for infleunce of winkler's foundation parameter on frequencey parameter A for 
guided-pinned beam 




pasternak foundation parameter 



Figure 5: Plot for Influence of Winkler Foundation Parameter on Frequency Parameter for 
Values of Pasternak Foundation Parameter for Various EC’s 

The influences of the foundation parametersy w andy p , on the stability parameter for guided-end conditions is 
shown in the Figure 5, The figures indicate that the stability parameter increases as the overall stiffness of the beam- 
foundation system increases. The overall stiffness of the beam- foundation system is an integrated resultant of the support 
stiffness, the foundation stiffness and the flexural rigidity of the beam. It is known that the flexural rigidity of the beam 
increases as X increases. It is obvious that the frequency parameter increases as the overall stiffness of the beam foundation 
system increases. 

The plots clearly show that while the Winkler foundation independently increases the frequency for any mode of 
vibration for constant values of warping and the Pasternak foundation parameters. Interestingly we can clearly observethat 
the effect of Pasternak foundation Parameter is to decrease the natural torsional frequency significantly for any mode of 
vibration and for constant values of warping and Winkler foundation parameter. 

The influence warping parameter Kon the stability parameter for guided-end conditions is shown in the Figure 6 
for all guided-end conditions. A close look at the results presented in Figure 6 clearly reveals that the effect of an increase 
in warping parameter K is to drastically decrease the fundamental frequency^. Furthermore, can be expected, the effect of 
elastic foundation is found to increase the frequency of vibration especially for the first few modes. However, this 
influence is seen to be quite negligible on the modes higher than the third. 
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plot for influence of warping parameter 'k' on frequency parameter X for guided- 
clamped BC for various values of Pasternak foundation parameter 




warping parameter 'k' 



plot for influence of warping parameter 'k' on frequency parameter X for guided- 
free BC for various values of Pasternak foundation parameter 




warping parameter 'k' 



plot for influence of warping parameter 'k' on frequency parameter X for guided- 
guided BC for various values of Pasternak foundation parameter 




warping parameter 'k' 



plot for influence of warping parameter 'k' on frequency parameter X for guided- 
pinned BC for various values of Pasternak foundation parameter 




warping parameter 'k' 



Figure 6: Plot for Influence of Warping Parameter on Frequency Parameter for Values of 
Pasternak Foundation Parameter and Winkler Foundation Parameter for Guided-End EC’s 
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CONCLUDINGS REMARKS 

In this paper, a dynamic stiffness matrix (DSM) approach has been developed for computing the natural torsion 
frequencies of long, doubly- symmetric thin-walled beams of open section resting on continuous Winkler-Pasternak type 
foundation. The approach presented in this paper is quite general and can be applied for treating beams with non-uniform 
cross-sections and also non-classical boundary conditions. Numerical results for natural frequencies for various values of 
warping and Winkler and Pasternak-foundation parameters are obtained and presented as graphical form showing their 
parametric influence clearly. From the results obtained following conclusions are drawn. 

• An attempt was made to validate the present formulation of the problem for guided-end conditions. There is very 
good agreement between the results, the general dynamic stiffness matrix defined by Eq. (33) and Eq. (47) of 
Yung-Hsiang Chen [15], for Euler-Bernoulli beam is observed to be same as Eq. (19) and Eq. (20) in this paper, 
but for only difference that the axial force was not included in the present paper. 

• The influences of the foundation parametersy w , y p and warping parameter K, on the stability parameter k for 
various supporting conditions are shown in the Figures (5 and 6). The second foundation parameter y p , tends to 
increase the fundamental frequency for the same Winkler constant]^, . The effect of y p , may be interpreted in the 
following way: A simply guided-pinned beam, which is the weakest as far as stability (k = 1.90) is concerned, 
acquires the stability of a beam which is guided at one end and free at other end (k = 5.22), by increasing the 
shear parameter of the foundation especially for the first mode. However, this influence is seen to be quite 
negligible on the modes higher than the first. Also, it is found that the effect of an increase in warping parameter 
K is to drastically decrease the stability parameter^. 

It can be finally concluded that for an appropriately designing the thin- walled beams of open cross sections resting 
on continuous elastic foundation, it is very much necessary to model the foundation appropriately considering the Winkler 
and Pasternak foundation stiffness values as their combined influence on the natural torsional frequency is quite significant 
and hence cannot be ignored. 

NOMENCLATURE 



Table 2 



T s 


St. Venant’s torsion 


T w 


warping torsion 


T 

1 n 


Total Non-Uniform Torsion 


K 


Modulus Of Subgrade Reaction 


P 


Pressure 




Shear Modulus 


0 


Angle Of Twist 


G 


Modulus Of Rigidity 


Cs 


Shear Constant 


M 


Twisting Moment In Each Flange 


h 


Distance Between The Center Lines Of The Flanges 


h 


Moment Of Inertia Of Flange About Its Strong Axis 


u 


Lateral Displacement Of The Flange Centerline 


C w 


Warping Constant 
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E 


Young’s Modulus 


p 


Mass Density Of The Material Of The Beam 


b 


Polar Moment Of Inertia 


K w 


Winkler Foundation Stiffness 


K p 


Pasternak Layer Stiffness 


z 


Distance Along The Length Of The Beam 


CO 


Torsional Natural Frequency 


K 


Non-Dimensional Warping Parameter 


Yp 


Non-Dimensional Pasternak Foundation Parameter 


Yw 


Non-Dimensional Winkler Foundation Parameter 


A 


Non-dimensional frequency parameter 


6»(z) 


Variation Of Angle Of Twist 0 
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